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Abstract 

In this paper we present a method for reducing the three dimensional Schrodinger equation to 
study confined metallic states, such as quantum well states, in a multilayer film geometry. While 
discussing some approximations that are employed when dealing with the three dimensionality 
of the problem, we derive a one dimensional equation suitable for studying such states using an 
envelope function approach. Some applications to the Cu/Co multilayer system with regard to 
spin tunneling/rotations and angle resolved photoemission experiments are discussed. 
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I. INTRODUCTION 



There has been much recent interest in magnetoelectronic devices due to their poten- 
tial applications as miniaturized computer memory components and high speed analog de- 
vices^. The ability to synthesize systems with artificial structures has grown remarkably 
over the past decade. For example, it has been shown recently that molecular beam epitaxy 
techniques allow for the growth of independent ferromagnetic layers on a semiconducting 
substrate, in cases such as Fe on InAs(lOO)^ and MnGa on GaAs(lOO)^. The search for new 
device materials and making optimum use of such new devices will greatly benefit from an 
accurate understanding of the underlying quantum mechanical processes involved in electron 
transport as the dimensions of the device approach the wavelength of transmitting electrons. 
Recent experimental studies of spin dependent, hot electron transmission such as those de- 

n ft 

scribed in Ref. U and resonant tunneling through two discrete states (Ref. |5|), have raised a 
number of interesting issues related to the ferromagnetic and insulating materials used, the 
nature of the electronic states that are involved in transmission and enhancements in spin 
filtering effects. Apart from the first principles based attempts which can be quite tedious, 
most theoretical studies of these spin dependent effects have used free electron band struc- 
tures and other simplifications in the metallic as well as in the insulating regions. Our work, 
though motivated by free electron approaches such as those introduced by Slonczewski^, is 
an attempt to bring out a more realistic lateral dependence of the electronic states under 
consideration. 

This paper is organized as follows. First, starting from the three dimensional Schrodinger 
equation, we proceed to derive an envelope function approach suitable for multilayered films. 
This procedure will go beyond the free electron methods that have been commonly used in 
the past, making use of more realistic wave functions, but avoiding a full pledged ab-initio 
calculation when studying such systems. We introduce an approach which incorporates the 
two dimensional Bloch wave vector k|| and show that the associated parallel band structure 
characteristic of the material being used, plays a major role in perpendicular transmission. 

Second, We will study spin tunneling and rotation effects (to be defined later) in a 
multilayer system with two ferromagnetic layers separated by a nonmagnetic metal such 
as Co/Cu/Co. We will also address some issues related to angle resolved photoemission 
and inverse photoemission experiments focused on confined states in metallic multilayers. 
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Although there have been attempts to explain such spectra using a Phase Analysis Model 
(PAM) with some success^, a better understanding of the multilayering, tunneling, spin 
rotation effects etc. is necessary, for example, to go beyond simple situations and study 
more complicated heterostructures of different sizes and shapes. For such systems, quantum 
mechanical calculations of spin-dependent transmission using a simple free-electron model 
have already been preformed^*- However, there is a growing need for more realistic, yet 
simple enough calculations in order to understand new and small structures that are being 
synthesized today. 

II. MODEL AND APPROXIMATIONS 

There have been some recent attempts to evaluate the effects of the two dimensional 
(planar), metallic periodic structures on the confined states in various devices using simpli- 
fied models 10 . Some of these have focused on a Fourier space description of the one particle 
Schrodinger wave function. We will be examining some of the assumptions made in these 
Fourier as well as real space models about the confined states, the effects of the planar 
regions on them and the importance of spin asymmetry on the spin filtering process. 

Our model consists of multilayered slabs of different materials sandwiched together to 
form a device. For example, the device may contain several layers of Cu sandwiched be- 
tween two ferromagnetic slabs of Co as in Co/Cu/Co. We will choose the z direction to be 
perpendicular to these slabs and the (x, y) to be parallel planes consisting of slabs. These 
can also be labeled longitudinal and transverse directions respectively. There are several 
simplifications that are usually made in attempts to calculate properties of such structures. 
When approximations are made to the wave function near an interface, it is a common 
practice to separate the transverse (x, y) dependence and the longitudinal (z) dependence 
sec Ref. Q). Let us carefully consider the conditions on the one particle potential and 
the tunneling states that lead to such a description of the problem at hand. An arbitrary 
eigenstate here can be expressed in the (planar) Bloch form as 




(1) 



G 



and the one particle potential U(x, y, z) can be expressed as 
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U{x,y,z) = Y^ V G n (z) exp{iG\\ ■ ry) (2) 

to elucidate the periodic nature of the potential and the Bloch-like form of the wave function 
in the planar (i.e. parallel) direction. Here Gm refers to a planar reciprocal lattice vector, 
while r|| = (x, y). 

It is possible to analyze the effects of the planar states on the perpendicular behavior in 
several different ways. First, using the Fourier coefficients defined above, the complete one 
particle solution may be expressed as, 



~^^ C Z^+ £ ^ lrG ; | W^ | W = {^-^(k|| + G||) 2 -^oW}^ l 1 l | W (3) 
G n^ G ii 

Note that unless the (parallel and z direction) potential coefficients V G ^_ G > (z) are weak com- 
pared to the relevant energy scale of the problem, the above equation couples the Fourier 
coefficients of the wave function through the potential coefficients and hence does not neces- 
sarily yield exponentially decaying solutions in the perpendicular (z) direction for the wave 
function even if the energy E satisfies the condition, 

£<|^(k|l + G||) 2 + ^) (4) 

for all values of z. 

Consider the wave function as defined in Eq. An assumption that is usually made 16 
when looking for such solutions is the following "separability" condition: 

* = £'(z,3/)0». (5) 

The assumption of the wave function separating into planar (£ ) and perpendicular (0 ) 
parts is equivalent to having Cq,Jz) = <j> (z)Dq,,,, where -Dg n Fourier coefficient has no z 
dependence. These ideas can be expressed in terms of the one particle potential U (x, y, z) 
and its Fourier expansion given by Eq. (j2J. It is clear that if this potential satisfies the 
(additivity) condition, U(x,y,z) = Ui(x,y) + U 2 (z) with no coupling between the planar 
x, y and perpendicular z dependencies, then the above separation of variables in the wave 
function can be justified. In such situations, the Hamiltonian H also becomes additive as 
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H(x, y, z) = Hi(x, y) + H 2 (z) and a simple quantum well equation, 

<P'(z) = Ec<I>'(z) (6) 

can be obtained for the (j)'(z). 

However, we argue that the above assumptions are too restrictive when one is dealing 
with states that have p, d or / character in local orbital angular momentum. The coupling 
of the three directions in the wave functions have to be dealt with more rigor, since such 
atomic wave functions are unlikely to satisfy the separability condition expressed in Eq. (j3J) . 
Significant corrections will be necessary if the separability condition is used as a starting 
point for better calculations. Hence we have sought a different starting point for carrying 
out the reduction of the three dimensional Schrodinger equation. The envelope function 
approach discussed below provides an ideal and formal platform for handling the situation 
at hand. 

III. ENVELOPE FUNCTIONS AND THE FULL PROBLEM 

With the advent of methods related to semiconductor quantum wells, the simplistic the- 
ory surrounding them has been so successful that sometimes it is easy to undermine its 
connections to the well understood regime of weak perturbations in bulk crystals. However, 
since the potentials employed in quantum wells are strongly perturbed at various boundaries, 
some formal justification seems necessary in order to use simple quantum well equations for 
multilayer systems. In fact, some applications of quantum well based techniques to semi- 
conductor heterostructures have been justified using an envelope function approach 14 . 

In principle, the relevant many particle Hamiltonian in all the different regions of the 
hetero junction carries all the necessary (interaction) information. When combined with the 
proper boundary conditions, its appropriate eigenstates can be used to describe, for ex- 
ample, tunneling in such a device. However, this problem is highly nontrivial and various 
approximations are sought in order to simplify it. First principles methods, such as those 
based on the density functional theory, can be utilized for this purpose but less complicated 
approaches that can reduce the computational burden are quite attractive. The envelope 
function method introduced by Bastard 12 is one such approach. Here the real space equa- 
tions satisfied by the envelope functions were equivalent to the k • p method of Kane 15 with 



rjf_^_ 

2m dz 2 



+ U 2 (z) 
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the band edges allowed to be functions of position. For hetero junctions with planar metallic 
regions, the applications of such ideas can be clarified and presented in a relatively straight- 
forward way starting from the nonrelativistic three-dimensional Schrodinger equation: 

V(x,y,z) = E*(x,y,z). (7) 

As is commonly done for (nonrelativistic) itinerant systems, the spin dependence in the 
Hamiltonian has been absorbed into a spin dependent potential U a , and the two equations 
representing the up and down spins can be solved separately. Hence, from now on, we 
will drop the spin index a and focus on the reduction of a single, one particle Schrodinger 
equation. In the parallel direction, the metallic as well as insulating regions are assumed to 
have perfect, two dimensional crystalline order, giving rise to extended electronic states with 
well defined parallel Bloch momenta frk||. We also assume perfect (parallel) lattice vector 
matches at various interfaces. 

An important point to note here is that in these problems that involve heteroj unctions, 
there are (at least) two relevant length scales; namely the interatomic length scale and the 
scale associated with the (confining) structures. The envelope function may vary on the 
latter length scale or on some combination of the two, which is to be determined later. 
Based on this argument, one can expect a nontrivial envelope function, when it exists, to 
modify the rapidly varying atomic wave functions. We regard this as our starting point and 
express the full problem (ignoring the spin dependence) as: 



-—V 2 + U a (x,y,z) 
2m 



2m dz 2 



+ 



h 2 



d 2 d 2 
+ 



2m \dx 2 dy 2 



U(x,y,z) 



£(x, y, z)<j>(z) = E£(x, y, z)<j>{z). 



The function £(x, y, z) can be thought of as a wave function with rapid variations on the 
atomic scale that has the two dimensional Bloch character, while <f)(z) is an envelope function 
as described above. An important point to note here is that we do not make an assumption 
on separability as in Eq. (jSJ). The existence of a nontrivial envelope function, as identified 
by the above equation, will be used as a prerqusite for the existence of quantum well or 
other confined states. Our search is for envelope functions, as defined above, that are likely 
to arise due to the confining structures. One can question the validity of such an expression, 
and similar forms have been suggested^ such as an expansion using products of Bloch states 
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and envelope functions. Here, we use the above form and associate an eigenvalue En with 
the function £(x,y, z) through the following eigenvalue problem. 



2m 



d 2 



+ 



d 2 d 2 
+ 



+ U (x, y, z) £(x, y, z) = E\\£(x, y, z) 



(9) 



dz 2 \ dx 2 dy 2 

We note that this is similar to how two dimensional band structure is calculated for thin 
films. Although for simplicity, we do not address interface roughness and other similar issues 
in this paper, the potential U (x, y, z) in Eq. (|SJ) can be modified to include such effects with 
appropriate changes in the boundary conditions (and a z dependent potential in Eq. JHJ). 
With the above definition of E\\ (subband energy), the full Schrodinger equation can be used 
to obtain a differential equation for the envelope function <f>(z) in the following manner. 



2 r 



h 
2m 



€(x,y,z] 







d 



dz 2 ^ z) + 2 Wz^ z) d-z^ z) 



{E-E\te(x,y,z)<i>{z) (10) 



The second term on the left hand side of the above equation contains a product of first 
order derivatives and describes some coupling of the lateral (i.e. planar) and longitudinal 
(i.e. perpendicular) coordinates, in addition to any coupling that is already contained in 
En. This coupling term can be simplified using averages of £(x,y,z) over the planar x, y 
coordinates. In general, this will result in a z dependent term and a second order differential 
equation for <p(z) as, 



with 



2m 



Pu(z) 



dz 2 



d_ 

dz 

M 



(E — E\\)(j){z) 



(2 J Jdxdy r^)/(/ jdxdy CO- 



(11) 



(12) 



The double integral over the planar coordinates x, y has to be carried out over a suitable 
2 dimensional unit cell. This is a mathematically rigorous result, based solely on the as- 
sumptions stated previously, and the theory at this level cannot and should not distinguish 
between metals, semiconductors or insulators. 

Now we can look for possible simplifications to Eq. ([lip by monitoring the properties 
of P\\(z). Note that if the Bloch function £(x, y, z) has an oscillatory z dependence, then 
P\\(z) will be purely imaginary. However, the imaginary part of P\\(z) (Im P{z)) is directly 
related to the (x, y) averaged flux, J(z) along the z direction. Hence, we can make use 
of (steady state) flux conservation which leads to dJ(z)/dz = 0; i.e. the conservation of 



flux implies that Im (Pm) has to be independent of z. When P\\{z) has a non negligible 
real part, we cannot make use of the above argument, and the type of confined states that 
we are searching for (through Eq. (|14|) ), will not exist. However, note that even for this 
situation, we have achieved the reduction of the 3 dimensional Schrodinger equation to a 
one dimensional one. Now setting P\\(z) = Q\\ (z-independent), the substitution 

<j){z) = exp{-zQ\\/2)C{z) (13) 

can be used to eliminate the first order derivative leading to a familiar equation, similar 
to Eq. ©. The function ((z) appearing here has a different interpretation, as a part of 
an envelope function, and the boundary conditions in the quantum well problem should be 
applied to the envelope function <p(z) or ((z). Note that that if ((z) and its derivative are 
continuous across various boundaries along the z direction, then similar properties can be 
established for d>(z). 



^ [£ CM 



dz 2 



h 2 Ql 

E-Eu{kn)- 



C(z) = EcC(z) (14) 



am 

The above equation illustrates several important points that are often overlooked or 
misinterpreted in simplistic quantum well and free electron approaches. As observed in Ref. 
licl ). a possible reduction in tunneling rates can be expected due to the lateral variation of 
the wave function. Although our method is somewhat different, note that Eq. (4) of Ref. 

n 

(JlfD has a (x, y) averaged term similar to Q in our discussion, which affects, for example, 
decay rates associated with s, p functions differently. The subband structure, i.e. £ l i|(k|i), in 
the multilayers affect the confined eigenstates and eigenvalues. The envelope functions are 
also affected by the averages of the parallel Bloch functions through Qn. For each 'confined 
energy level' with energy Eq, there exists a continuous subband of states that share the 
same perpendicular wave function <j>{z) but differ in parallel Bloch momentum kjj, En and 
£(x,y,z). Note that E\\ can depend on the thickness of a given multilayer film and carries 
information not only about the planar structures, but also about the longitudinal coupling, 
following our definition through Eq. (jUJ). Finally, the energy E of a given electron in such a 
quantum mechanical state depends on all of the above. 
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IV. APPLICATIONS - CONFINED STATES IN METALLIC MULTILAYERS 



We are now ready to apply our model to test devices. The first device consists of sev- 
eral layers of (nonmagnetic) Cu, sandwiched between two slabs of (ferromagnetic) Co (i.e. 
Co/Cu/Co) which will be used to discuss spin transmission/rotation effects. For an incident 
electron of energy E tota i, we obtain the following electronic perpendicular momenta, 

hk u = pm{E total -E\\)-h 2 QyA (15) 
hkf z = pm{E total -E\\-/\)-h 2 Qy± (16) 

where fcf (Ar^) is the perpendicular wave vector for majority (minority) electrons in the ferro- 
magnetic, metallic regions where the lateral effects have been taken in to account using the 
ideas developed in the previous sections. Here A is the spin splitting in the two dimensional 
bands assumed to be kn independent. 

In this device, we can rotate the magnetization of the right (R) ferromagnet by an angle 9 
with respect to the magnetization of the left (L) ferromagnet. The spin rotation is introduced 
at the boundary z = inside the spacer (Cu) layer for simplicity. For this device, solutions 
to Eq. (|14jl in various spatial regions (L and R, subdivided into a, c, d, e) take the form, 



SaJ 
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(17) 



Finally, we impose the boundary conditions at z±z 2 (where 2z 2 is the thickness of the center 
slab) and at z = for spin rotations. 







C c L T (-^ 2 ) 


.CaV^2)_ 
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Cf T (o) 
CS(o) 



C L T (o) 



c5(o) 



/ 



cos(#/2) sm(0/2) 



(18) 



is the 



and identically for the first derivative. The matrix S{6) = 

\ -sin(6/2) cos(9/2) 

spinor transformation, where 8 7^ is tied to the spin rotation effects as discussed below. 
We can then fully determine the transmission coefficients utilizing numerical techniques. 



A. Spin Transmission and Rotations 

An interesting application of the above ideas has to do with ballistic transport and spin 
rotation effects. When polarized electrons are transmitted from one region to another with 
a different polarization axis, they experience a spin-torque and a transfer of angular mo- 
mentum to the new medium-*- These ideas have now been demonstrated experimentally, 
for example, through the phenomenon of giant magnetoresistance where large current den- 
sities flowing perpendicular to the films have been observed in reversals of magnetization^, 
and spin precession^. This field is an emerging one related to 'magnetoelectronics' and 
many new experiments are expected to be conducted on spin transmission/rotation effects 
in magnetic multilayer systems. 

We do not wish to focus on the mechanisms of spin transfer but simply use spin rotation 
angle as an input to our calculations and obtain the corresponding transmission coefficients in 
a device consisting of two magnetic films separated by a nonmagnetic one (i.e. Co/Cu/Co). 
Spin rotation effect is introduced, by hand, at the center of the nonmagnetic film (Cu) 
for simplicity as has been done previously^. However, unlike in the previous studies, the 
underlying band structure and lateral effects have been taken into account by using the 
theory of confined states developed here. These band dispersions play a crucial role in 
determining the spin dependent transport properties. For example, when such confined 
states in a ferromagnetic film are located in a gap of minority spin states, the system can 
act as an almost perfect spin filter . The k-dependent transmission coefficients, and T^, 
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k„ (a; 1 ) 

FIG. 1: Transmission coefficients T^,T^ as a function of kn along TX in the second device 
(Co/Cu/Co) for a relative spin orientation 9 = j. At the zone boundary X, the device filters 
out the 'minority' spins. Away from the zone boundary both 'majority' and 'minority' spins are 
transmitted. 

that were introduced in the previous section can be used through a Landauer type formula^, 

G = jY,nh) (19) 

k ll 

to obtain the conductivity in the quantum well problem under discussion. Here we calculate 
transmission coefficients along TX using selected subbands, as defined in Eq. f!20|) . 

#ll(k||) = E§ + W{1- cos(k x a)}. (20) 

in the transmission device (Co/Cu/Co) for the following set of parameters (Fig. for 
illustrational purposes: W(Co)= 0.07, W(Cu)= 0.06, A = 0.19 (all in Rydbergs), z 2 = 13.6, 
a= 6.8 (all in Bohrs). 

At a given (total) energy of the incoming electron, and a given relative spin orientation of 
the Co films, transmission coefficients T^,T^ have been obtained from Eqs. ()17|) and (|18|) 
as a function of the parallel Bloch momentum, kn. From this figure, it is easy to see the 
effects of different types of Bloch states on the tunneling. Near the zone boundary X, for 
the given relative spin orientation 6, we see the device filtering out 'minority' spins, while at 
other k-points along TX, a mixture of both 'majority' and 'minority' spins are transmitted. 
This is a direct result of the upward dispersion of the selected subbands along FX and the 
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spin splitting in Co, pushing the minority band closer to the given (total) energy of the 
electron. 

B. Energy Spectrum and Photoemission 

Angle Resolved Photoemission (ARP) techniques have become a useful tool for studying 
electronic states and their (band) dispersions in structures with two dimensional (planar) 
Bloch symmetry. In ARP experiments photons are used to eject electrons from occupied 
states while in inverse ARP, photons are ejected, when the above photoemission process 
is reversed. Photoemission is a many body phenomenon and an evaluation of the spec- 
tral function and appropriate self-energies yield the full photoemission spectrum and peak 
widths. However, here we follow a simpler approach and focus on one particle energy states 
and the changes introduced by the multilayering. In this approach, the energy E that is 
measured for the photoemitted electron can be associated with the full Schrodinger equa- 
tion (j2J). Existence of an envelope function as defined in Eq. (jHJ and satisfying the imposed 
boundary conditions are necessary for confined states. The energy En as defined through 
Eq. © has to do with the ordinary subband states due to the periodic potential in the 
planar directions. If the multilayering does not play any role for a given energy E at a given 
k||, then En and E have to be identical when these solutions exist in Eqs. (J7j) and (0). In 
such cases, <p(z) turns out to be a simple multiplicative constant and Eq. (fTUj) is consistent 
with this scenario. However, when confined states exist, they may alter the usual dispersions 
En that are observed in their absence. According to the theory developed above, confined 
states can be identified as states for which solutions to the differential equation (j!4|) exist 
for a given energy E of the electron. It is also important to realize that, for a given energy 
E, these states may not exist for all values of k|| along a given direction of the 2 dimensional 
Brillouin zone. 

In general, the dispersion of the total energy of the electron E, i.e. the energy of the 
bound electron as a function of k||, are found whenever confined states exist with appropriate 
as solutions to Eq. (JHJ. In the present context, with boundary conditions appropriate 
to stationary states (determined by the film geometry), we obtain a set of discrete energies 
(E n ) and states for the one dimensional quantum well or barrier problem where 
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E=E n + E ]] (k ]] ) + -^. (21) 

Hence, the observed energies E in an ARP experiment for confined states will depend on 
the existence of these discrete E n s and the corresponding ( functions satisfying the relevant 
boundary conditions. As in simple quantum well problems, solutions of odd and even parity 
may sometimes be associated with E n and these will depend on the film geometry, interfaces 
and growth conditions. We plan to address these as well as effects due to film thickness 



(observed in Refs. 



21 



in a future paper. 



V. CONCLUSIONS 



Using an envelope function approach, we have reduced the three dimensional Schrodinger 
equation to a one dimensional one and analyzed the lateral effects of metallic multilayer 
films through a hypothetical subband structure that varies with the two dimensional Bloch 
vector. Ballistic transmission in a magnetoelectronic device consisting of Co/Cu/Co has 
been examined and shows a clear dependence on the Bloch vector k||. Energy spectrum 
of confined multilayer states that may be observed in angle resolved photoemission studies 
are discussed. Also, this method of reduction of the Schrodinger equation with appropriate 
modifications may be used in other confined geometries such as nanoparticles and wires. 
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